Self-organized Criticality and Absorbing States: Lessons from the Ising Model 
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We investigate a suggested path to self-organized criticality. Originally, this path was devised to 
"generate criticality" in systems displaying an absorbing-state phase transition, but closer examina- 
tion of the mechanism reveals that it can be used for any continuous phase transition. We used the 
Ising model as well as the Manna model to demonstrate how the finite-size scaling exponents depend 
on the tuning of driving and dissipation rates with system size. Our findings limit the explanatory 
power of the mechanism to non-universal critical behavior. 
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Self-organized criticality (SOC) refers to the sponta- 
neous emergence of critical behavior in slowly driven dis- 
sipative systems 0, |3l ■ Most models are defined on lat- 
tices with local particle numbers Zi and thresholds zf. 
They are driven discretely in time by increasing Zi at 
randomly chosen positions i until such an increase leads 
to Zi > zf somewhere in the system. Particles then top- 
ple to neighboring sites and can trigger avalanches of 
local redistribution propagating through the entire lat- 
tice. Dissipation typically takes place at the boundaries, 
where particles leave the system. When an avalanche 
has finished the model is driven again 0. The resulting 
avalanche size distributions obey simple scaling. 

In models displaying absorbing state (AS) phase tran- 
sitions :3-] a tuning parameter, such as the overall particle 
density, controls a transition between an inactive phase 
and a phase where activity in the system continues in- 
definitely. 

From the first introduction of SOC in 1987 it 
was believed that SOC models manoeuver themselves to 
the critical density between similar inactive and active 
phases. Tang and Bak suggested in 1988 that the density 
of "lattice sites on which z > Zc [■■■] may be viewed as the 
order parameter for this critical phenomenon" 4]. Such 
an identification of the activity with the order parameter 
implies a link to absorbing state phase transitions. 

This link was formalized and made explicit about 10 
years later d, 0, ■ Dickman et al. introduced pe- 
riodic boundaries to SOC systems, thereby turning them 
into AS models. Measuring the exponents characterizing 
the spreading of perturbations |3, 13 or the roughness 



of the associated interface models has been 

observed that at the critical density the closed-model be- 
havior resembles that of open SOC models 

The resulting interpretation of SOC is obvious 
Activity eventually leads to dissipation at the bound- 
aries, which in turn reduces the particle density to below 
the critical value. Driving takes place whenever quies- 



cence has been reached. SOC models therefore hover 
around the critical point, being pushed forth into the ac- 
tive state by driving and pushed back into the quiescent 
state by dissipation. 

With this simple picture in mind one arrives at an 
equation of motion for the particle density in the system 



eqtiatioi 



ds 



C(s) = h- Pa{s)t , 



(1) 



where s is the time, h is the driving rate and e is called 
the (bulk) dissipation rate. The activity pa is the order 
parameter, defined as the density of active sites, Zi > zf, 
in the active phase. We will refer to this interpretation 
of SOC as "the AS approach" . 

Clearly, the driving h must be very slow compared to 
the dissipation pae. Otherwise particles would be added 
while the system is active, leading to a fiuctuating activ- 
ity rather than distinct avalanches. The proponents of 
the AS approach point out that h, e and h/e have to be 
tuned to zero in order to achieve the desired separation 
of timescales 0|. While the definitions of SOC models 
typically restrict dissipation to boundary sites and re- 
sult in diverging avalanche sizes in the thermodynamic 
limit, leading to appropriately vanishing e{L) and h{L), 
so far no statement has been made as to how the limit- 
ing behaviour is approached. But this turns out to be the 
all-important piece of information: The finite-size scal- 
ing (FSS) behavior, the only scahng available in SOC, 
depends entirely on the scaling of the driving and dissi- 
pation rates with system size. Choosing h{L) and e{L) 
freely, arbitrary scaling behavior is produced. 

In the following the relation between the scaling of h 
and e and the resulting FSS is analyzed, using the two 
dimensional Ising model as an example. However, the 
analysis is generally applicable and works equally well for 
standard SOC models and their AS counterparts, which 
is confirmed by simulations of the Manna model |llL Il5l | . 
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Translating Eq. iQJ into magnetic language, ( corre- 
sponds to the inverse temperature /3 and the activity pa 
to the modulus of the magnetization density |r7i|. The 
parameters h and e become cooling and heating rates, so 
that the temperature T is increased for large magnetiza- 
tions and lowered otherwise, 

^^(3is)^h-\mis)\e. (2) 

The resulting model is an Ising model where the temper- 
ature is dynamically adapted according to the equation 
of motion Therefore, the configurations are not sam- 
pled with Boltzmann-weight and the resulting "dynami- 
cal ensemble" is not canonical. However, by multiplying 
(O by a small pre- factor, corresponding to rescaling the 
time, the distribution of temperatures can be made ar- 
bitrarily narrow. For the sake of the following analysis, 
it is assumed that this "dynamical Ising model" is well 
characterized by a single effective (reduced) temperature, 

teS- 

For the FSS analysis presented below we choose the 
approach of /i, e ^ to leading order as 

h = hoL""^ and e ^ eoL"'", (3) 

where uj,k > 0. In the stationary state, (j^P) = 0, 
@ yields (|m|) = {ho/eQ)L'^-'^ with () denoting the 
average over the dynamical ensemble introduced above. 
Clearly one must choose uj > k. To attain the prescribed 
(|m|) (L) the system settles at the effective (reduced) 
temperature {T -T^) /T^ = tas{L) oc L"!/'' to leading 
order, see Fig. Via tes{L) all thermodynamic quanti- 
ties depend only on L, which can be mistaken for stan- 
dard FSS at temperature T = Tc. For the study of SOC 
models it is vital to understand the difference because 
SOC systems are always critical, wherefore FSS is the 
only scaling available. 

Around the critical point of a continuous phase transi- 
tion, the singular part of the free energy leads to a simple 
scaling behavior of the magnetization density jl^]. 

(|m|) = ^kf^L-f'/^VihtL'^n, (4) 

where t is the reduced temperature, negative in the low 
temperature phase (LTP) and positive in the high tem- 
perature phase (HTP). and kt are metric factors, and 
Y'{x) is a universal scaling function, which becomes de- 
pendent on the boundary conditions and the geometry of 
the system in the limit of small arguments, case (|5b|l be- 
low. There are three qualitatively different (asymptotic) 
regimes 
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where in the Ising model ^ — vd ~ 2/3. 
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FIG. 1: Finite size behaviour for the magnetisation \m\ versus 
the temperature T in the regular, 2D-Ising model for L\ — 80 
(squares) and L2 = 160 (circles). The dashed line shows 
Yang's solution p^. In the dynamical Ising model, the cool- 
ing and heating rates, h{L) and e(L) prescribe the magneti- 
sation |m(I/)| — h{L)/e(L) oc L'^~^ , indicated by horizontal 
lines. The system is forced to move to an effective tempera- 
ture te«{L), indicated by arrows. 

The first line describes the asymptotic behavior of the 
magnetization in the LTP, the second line represents FSS, 
and the third line describes the HTP. 

Setting (|m|) oc L'^^'^ , the different regimes of ||SJ) are 
accessed by three qualitatively different choices oi k ~~ lu, 
that is speeds at which (|m|) approaches zero: 

1) K — Lj > —l3/v ("too slow"): In this case the mag- 
netization approaches slower than in a standard 
Ising model kept at temperature T = Tc as the 
system size increases, so that Y' {kttcfi{L)L^/'^) oc 
j^K~ui+l3/v jg divergent in L. The only way to ob- 
tain a divergent Y'{x) is via (|5a|l . which requires a 
negatively divergent argument a; ^ — 00. The ef- 
fective temperature is therefore negative and scales 
like |ioff(i)|^ oc L'^"". Using ieff(i) oc L"!/^ leads 
to 

^JL^ (3/{uj- k)> V . (6) 

This implies that tcs{L) finally leaves the FSS re- 
gion, whose width scales like L"^/", toward the 
LTP. 

2) K ~ UJ ~ —P/v ("correct"): In this case Y'{x) re- 
mains constant, so that its argument either remains 
constant or vanishes, according to (|5 b|) . Thus 
tcs{L) decays at least as fast as L"^^", i.e. fi < v. 
To the order considered here the equality applies. 

3) K — u}< —jS/v ("too fast"): Y'{x) vanishes, follow- 
ing (|ScJ|. For Y'{ktt^s{L)L^/'') oc L^'-^^+P'" and 
tcs OC L^^/^ this implies 

7^ 
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provided that the denominator of Eq. lO is posi- 
tive. Hence the model leaves the FSS region toward 
the HTP. The special case of negative /Lt, implying 
divergent effective temperature, will be ignored. 

Crucially, only for k—uj — —[3/u (case 2)) does the model 
remain in the FSS region. To achieve this, h/e must be 
tuned exactly in the way the order parameter scales in 
a system displaying standard FSS, oc L'^^" while 

fixed at the critical temperature. In all other cases the 
scaling of the effective temperature eventually drives the 
model out of the FSS region: f/L vanishes in the ther- 
modynamic limit. Nevertheless, (T) converges to Tc, so 
that the correlation length 

C oc L"^/^ (8) 

diverges. With this scaling of ^ all observables will show 
standard finite-size scaling with v replaced by fi |29j . 

To illustrate the above analysis we performed simu- 
lations of an Ising model with dynamics as described: 
Using Metropolis updating, the absolute magnetization 
density is calculated after each scan over the lattice. Ac- 
cording to (121) a new temperature is then calculated to be 
used in the next sweep, T = —h -I- (|m|) e. Starting from 
T = 2.27, systems of size L = 40, 80. ..640 were updated 
at least 10^ times as transient and at least another 10^ 
times for statistics. 

Our numerical simulations fully confirm the above 
analysis: We observe the standard FSS exponents with 
V replaced by ^ for any reasonable choice of k — uj. The 
new scaling exponent /x (and Tc) can be determined from 
(T) - Tc oc L-^/t". Using it in an FSS analysis allows 
us to identify all standard critical exponents. Even with- 
out the knowledge of n three measurements, say a//i, 
13/ ^ and 7//^, are sufhcient to determine all exponents 
using standard scaling laws. This seems to defy com- 
mon sense, since the exponents are determined without 
referring to Tc, making it seemingly very attractive for 
investigations of continuous phase transitions. However, 
standard methods, such as FSS or analysis of the critical 
behavior at ^ <C T, are much more reliable and efficient: 
Not only is the above identification of t^s with the re- 
duced temperature questionable. More importantly, al- 
most any choice of the scaling of h and e leads to vanish- 
ing £,/L. One therefore simulates effectively independent 
patches of a lattice in a way that FSS effects remain im- 
portant. 

For two reasons the method is very sensitive to the 
choice of eo and in Eq. Q [l^l' Firstly, the ampli- 
tudes of the fluctuations in the effective temperature de- 
pend directly on h and e; choosing /iq and eq too large, 
the system destabilizes. One can estimate these fluctua- 
tions by analyzing (j^J and derive a lower bound for 
Secondly, if h and e are too small and initially place the 
system close to Tc, the scaling function reaches its asymp- 
totic behavior (generally ifScjl or (|5a|l 'l only for very large 
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FIG. 2: The scaling of the effective temperature tctt (filled 
symbols, full lines) and the susceptibility x (open symbols, 
dashed lines) in the Ising model for different choices of the 
exponent k and u = 1. The symbols are numerical simula- 
tions, the lines show the slopes expected from theory. 

system sizes. 

Fig. El shows the scaling of the effective temperature 
for the three qualitatively different choices of the driving 
exponents discussed above. The values of /j, and 7 derived 
from these data confirm the calculations. Depending on 
the choice of k — w, the value of immediately determines 
either /3, Eq. ©, or 7, Eq. 0. The FSS of specific heat 
and susceptibility produces the expected values of 
and 7//Lt. 

Since our interest in the AS approach is due to its 
proposed role as an explanation for SOC, we repeated 
the analysis for a variant of the one-dimensional Abelian 
Manna model This sandpile-like model has been 

used to exemplify the link between SOC and AS 
It is driven in the bulk and implements bulk dissipation 
[l8| . as suggested b y ffl . While the key equation Q 
has been confirmed UJl, the Manna model is neither as 
well understood nor as well-behaved as the Ising model. 
Unlike in the Ising model, the Manna model can get 
stuck when hitting an inactive state. This complicates 
the analysis especially for the third case discussed above. 
Nonetheless the effective particle density {() is certainly 
a function of h and e and therefore depends on their scal- 
ing with all the consequences laid out above. Indeed, the 
numerics agrees fairly well with the theoretical predic- 
tions. Most crucially, the activity as well as the scaling 
of the avalanche size distribution show a clear, immediate 
dependence on the choice of the two exponents u and k. 

One should note that introducing bulk-drive and - 
dissipation to SOC models does not result in a full cor- 
respondence to AS models. Firstly, there are important 
observables in SOC lacking a counterpart in AS and vice 
versa. For example, there is no obvious definition of the 
avalanche size in the active phase of an AS model. Simi- 
larly, the definition of the SOC-activity is somewhat arbi- 
trary in Abelian models ,19 .] , due to the lack of a unique 
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microscopic timescale, which is needed when taking tem- 
poral averages or measuring rates. Secondly, for AS mod- 
els Eq. (0J) contains the asymptotic conditional activity, 
while in bulk-driven SOC models, such as the Manna 
model presented above, the instantaneous activity enters 
into the equation of motion (Q) . 

The present study shows that the proposed explana- 
tion of SOC as "self-organized" AS criticahty [1,0, 13 im- 
plies non-universal scaling behavior dependent on k — w. 
Universal, that is dissipation- and driving-independent, 
scaling behavior cannot be achieved with the AS ap- 
proach. The question whether SOC models have univer- 
sal features is very important. Universality is a main jus- 
tification for studying simple models and for disregarding 
the details of the physical processes they describe. 

Despite the importance of this issue, it is still unclear 
whether SOC systems can be grouped into universal- 
ity classes; in fact, exponents can change due to small 
changes in the update rules (e.g. 113,01), SOC is 
notorious for its wide variety of critical exponents. Ac- 
cepting the AS approach this would be a consequence of 
implicitly setting the scaling of external drive and dissi- 
pation by the dynamical rules of the different models. 

However, there is also strong evidence in favor of uni- 
versality in SOC. Many changes of the detaile d dv namics 
do not affect the critical behavior [H El Q 113 • More- 
over, the ratio £,/L appears to remain constant in direct 
measurements of some models so that the "correct" 
FSS exponents are observed, which is in stark contrast 
to ip. 

At first sight, the observation of the same exponents in 
SOC and AS models (such as 0,0,0]) seems to support 
the case for the AS approach. But our analysis shows 
that the opposite is true: If the AS approach was deter- 
mining the behaviour of SOC models, it would almost 
certainly (apart from case 2)) produce exponents which 
differ from those observed in their AS counterparts. Bar- 
ring coincidence, observations supporting universality in 
SOC must therefore be taken as strong evidence against 
the AS approach explaining the critical behavior of SOC 
models. 

We have calculated the FSS behaviour of a system 
approaching its critical point through a feedback mech- 
anism between order parameter and tuning parameter. 
While scale-free distributions of responses such as those 
observed in the case of rainfall or earthquakes [28j |. 
can be produced by such a process, it only yields critical 
behavior strongly dependent on the detailed dynamical 
rules of SOC models. There would be no universality and 
robustness against small changes in the dynamical rules. 
While the AS mechanism in its present form may pro- 
duce further insight into potentially non-universal criti- 
cal phenomena as observed in field experiments, it fails 
to explain the apparent universality of SOC models. 
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